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In this work, we introduce an extra singlet pseudoscalar into the Type-III two Higgs doublet model 
(2HDM) which is supposed to solve a series of problems in the modern particle-cosmology. With 
existence of a light pseudoscalar, the h > jit excess measured at CMS and as well as the (g — 2), 
anomaly could be simultaneously explained within certain parameter spaces that can also tolerate 
the data on the flavor-violating processes T — py and Higgs decay gained at LHC. Within the 
same parameter spaces, the DM relic abundance is well accounted. Moreover, the recently observed 
Galactic Center gamma ray excess(GCE) is proposed to realize through dark matter(DM) pair 
annihilations, and in this work, the scenario of the annihilation being mediated by the pseudoscalar 


is also addressed. 
I. INTRODUCTION 


Even though in the framework of a minimal ex- 
tended Standard model (SM) with non-zero neutrino 
masses the leptonic flavor violation (LFV) process is 
negligible for the smallness of neutrino masses which 
are experimentally confirmed [i]. Therefore, a di- 
rect search for the LFV processes would provide an 
ideal probe for new physics beyond SM, or in other 
words any observational anomaly may hint us its ex- 
istence. Besides the B-factory, since its high energy 
and luminosity, LHC is definitely the machinery for 
the exploration. A search for LFV has been per- 
formed by the CMS collaboration via two channels 
h — pte and h + utp, and a 2.40 excess of the 
branching fraction BR(h > pr) = (0.8473'3°)%) 2] 
is reported. If one could account such an excess as 
an anomaly, there should be some mechanisms which 
are obviously beyond SM, to be responsible for it. 
The Type-III two Higgs doublet model (2HDM) is 
one of them, because in the model a flavor-violating 
Yukawa interaction exists which may contribute to 
the LFV at tree level. The model has been explored 
extensively to study this phenomenological 
observation. Furthermore, the Yukawa interaction 
contributes to the muon g — 2 via one loop dia- 
grams and thus would provide a possibility to ex- 
plain the (g — 2), discrepancy [12]. Mean- 
while, the flavor-changing Yukawa interaction would 
induce a substantial contribution to the radiative de- 
cay T — uy, thus the flavor-changing Yukawa inter- 
action might be rigorously constrained by the avail- 
able experimental data [9] [TT] [12]. 


One of the main characteristics of the 21st cen- 


tury is that the cosmology has already become an 
accurate science and the corresponding observation 
must be combined with the precise measurements 
and new discoveries at the facilities on the Earth to 
testify the standing theories. The identity of dark 
matter (DM) and the interaction which determines 
the behavior of DM particles are the key point and 
searching for them is the most challenging job for 
both experimentalists and theorists of high energy 
physics and cosmologists. 

Recently, the Fermi Large Area Telescope data 
show an excess of gamma-ray at energy of a few GeV 
coming from Galactic Center (GCE) [321]. To ex- 
plain the observation, it is suggested that annihila- 
tion of DM particles weighing 30— 70 GeV into bb is 
responsible for the GCE [22}{24]. Even though there 
exist other proposals to explain the excess, such as 
a population of millisecond pulsars (MSP) [25} [26] 
which might be responsible for GCE, it is not easy 
to explain the energy spectrum and spatial distri- 
bution of the GCE [27] [28]. Thus in this work we 
discard the astrophysical source explanation 30) 
and focus on the DM scenario. 

The dwarf galaxies are considered to be the clean- 
est sources for detecting gamma rays produced by 
DM annihilations, thus the data on gamma ray ob- 
served at Reticulum II [31}{33] may imply existence 
of abundant DM at our galaxy. To be sure, we need 
to compute the cross section of DM annihilation in a 
model of particle physics. Meanwhile other cosmo- 
logical phenomena must also be concerned, namely 
the DM annihilation cross section required by the 
new data should be of the same order as that deter- 
mined by the thermal DM relic abundance. 

However, the original Type-IIJ 2HDM does not 


provide a natural explanation to the DM particle 
annihilation cross section. Thus we need to extend 
the model which can also accommodate the DM 
annihilation. It is note that a pseudoscalar could 
mediate the annihilation of dark matter(DM) pair, 
meanwhile due to the small momentum transfer at 
t-channel, the interaction between the DM particles 
coming from the outer space and the nuclei in the 
detector is not affected by the existence of the new 
pseudoscalar, so that the DM particles may evade 
the direct search at not-much sensitive detectors. 
This idea has been implemented in various mod- 
els [24 [34}{48]. It further motivates us to consider 
the DM explanation to the GCE, and introduce a 
Dirac fermion field serving as the DM candidate. In 
this work, we introduce a pseudoscalar ap into the 
Type-III 2HDM. The pseudoscalar does not directly 
couple to the SM particles, but it slightly mixes with 
the CP-odd Higgs which exists in the original Type- 
IHI 2HDM, thus it would effectively couple to SM via 
this mixing. Therefore the pseudoscalar can mediate 
an effective interaction between the DM yx and the 
SM fermions bb. 


Moreover, its introduction may bring up two more 
advantages as following: 


1. An extra pseudoscalar would open a new decay 
channel for the Higgs and thus affect the h > uT 
excess. 


2. There could be a mixing of the newly intro- 
duced pseudoscalar with the CP-odd scalar Ap in 
the original Type-III 2HDM, and hence induces a 
new contribution to the value of (g — 2),,. With in- 
crease of the mixing, the contribution of the new 
pseudoscalar would cancel that of the Heavy Higgs, 
thus in this extended model the theoretical predic- 
tion on (g—2),, can be decreased to a tolerable level. 


In this work the pseudoscalar also plays a role to 
explain the h — uT excess and the discrepancy be- 
tween theoretical prediction and data of (g — 2),. 


The paper is organized as follows: after this intro- 
duction, we discuss the new scenario where a light 
pseudoscalar is introduced to extend the Type II 
2HDM in section. |I| In the section. we investi- 
gate the relevant topics including the h —> uT excess 
observed at CMS, muon (g—2),, anomaly, the galac- 
tic center gamma ray excess (GCE) and the con- 
straints from the T —> py process and LHC Higgs 
data. The numerical results are presented in sec- 
tion. then the last section is devoted to our con- 
clusion and discussion. 


I. THE MODEL 


In this work, a Dirac fermion (x) of mass my 
stands as the DM candidate and a gauge singlet 
pseudoscalar ag is introduced to extend the Type- 
II two Higgs doublet model, where ag mediates the 
coupling between the dark matter and the SM par- 
ticles. The interaction Lagrangian reads 


Laark = —Yyyaoxiy?x ‘ (1) 


The pseudoscalar a9 mixes with the pseudoscalar in 
the original 2HDM and then couples to the SM par- 
ticles through the potential, given as[48] 


1 Aa 
V= V2HDM F 3a ag F Th F Vportal; (2) 
and 
Voortal = —iBao HÌ He + h.c. , (3) 


where B is a parameter of mass dimension, and the 
Higgs potential is{49] 


Vəupm = m1 HÌ Hy + oH Ho + (u3H! Ha + h.c.) 
Ny (mm) 154 (mm) +3 (HÌ) 
x (11) 12) = 7 (Him) (1m ) 
+ [(AsH] He + AsH] My + ArH} Eh) 
x (Him) Es hc. (4) 


We can explicitly rewrite Hı and Hə in the Higgs 
basis as 


Gt Ht 
A, = v+¢1+iG? , H= ( ob2tiAo ) ’ (5) 
V2 V2 


where Gt and G° are the Nambu-Goldstone bosons 
and H* and Ag are a charged Higgs boson and a 
CP-odd Higgs boson, respectively. Without losing 
generality, let us concentrate on the CP-conserving 
case, where ag does not develop a vacuum expecta- 
tion value (VEV) and all A; and u3 are set to be real. 
Then the potential is minimized to (H1) = v/V2, 
(H2) = 0, with v = 246 GeV. 
In the basis of (41,2), the mass matrix for the 

CP-even Higgs is M? whose elements are, 

Mĉ i = 2r1v7, 

dav? 
Miron = Migs + “5 + sv”, (6) 
Miz = Mig, = 60”. 


Diagonalizing the matrix, one obtains the physical 
CP-even states h and H (Mma < my) as 


go, \__[/ cosa sina H (7) 
ġ2 J \ —sina cosa h J?’ 
2Mi12 

Mion ~ Mhi 


with eigen-masses being 


tan 2a = 


1 
Ma H =3 Min + Minos 


+ JRun = Mea +4 Mia , 


and we consider the eigenstates h and H as the SM- 
like and heavy Higgs bosons respectively. 

The CP-odd Higgs Ag mixes with ag due to the 
potential Vyortar (Eq. B), and the mass matrix in the 
(Ao, ao) basis is 


m, Bv 
me e N (8) 
ao 
where m?o = igre + r4v2/2 — Asv?. Thus, the 


relation between Ap, ao and mass eigenstates A and 


a is 
(#)=( 35 aN (9) 
ao sin cos@é a 
with the mixing angle 
2Bv 


2 — 2 
MAg Mao 


6 = = tan™!( 


), (10) 


and the masses squares are 


1 
mZ A = [rs + m2 + Vm, — m2) + 4B2v2 


The parameter B in terms of ma,a and 0 can be 
expressed as 

B= zs (m% — m?) sin 20. (11) 

2v 

The effective coupling of the CP-even Higgs bosons 
to the SM W is igmwCpwwg"” with Chww = 
sina, Cyww = cosa, and the CP-odd Higgs bosons 
A(a) do not couple to W, i.e., Cajayww = 0. And 
the Vportal is recast in terms of mass eigenstates and 
mixing angle as 


1 
Vportal = ae (m3 = m2) [S40 aA+ oe (A? = a’) | 
x (sina h + cosa H). (12) 


The effective coupling of DM fermions to the medi- 
ator given in Eq. is simply expressed as, 


Laark = —Yy (cos 0 a + sin 0 A) xig? x. (13) 


The Yukawa interactions in the extended Type-III 
2HDM are 


LYukawa = -Qi Vå uiyidR = O Vea Hp dh 
— Qi iyluh — Qi Ap uh 
— Li, Hiyter — Li, Hopi eh, (14) 


where Q = (ut, Voxmdz)’, L= (VMNSYL, eL)” and 
Hi, stands for iogH*. Voxm(Vuns) is the Cabbibo- 
Kobayashi-Maskawa (Maki-Nakagawa-Sakata) ma- 
trix. The general 3-by-3 complex matrices p? induce 
the Higgs-mediated Flavor Changing Neutral Cur- 
rent (FCNC). In the mass eigen-basis of the Higgs 
bosons, the Yukawa interactions are recast as 


LYukawa = —Ygij frit fri om Dri(Vinspe) H er; 
—ti(VexmpaPr — pl Vorm Pr)” H* d; 
+h.c., (15) 


with ọ = h, H, A,a, f =u, d, e,v, and 


ij 


mi p 
fa f 
Yhij = — sin ad;; + V2 COS Q, 
mi pi 
YHij = a COS Adj; — A sina, 
. ij 
- “E cos 0, (f = u), 
YAij = ip 
z 00s, (f =d, e), 
ip 
— sing, (f = u), 
Yaij = E (16) 
Ea sind, (f =d, e), 


where the couplings y4;; and Yaij exist in the new 
Feynman rules and are accompanied by y5. For in- 
vestigating (g — 2)„ excess, we do not need to in- 
voke the so-called Cheng-Sher ansatz for p? [50] 
since the corresponding parameter space is highly 
restricted [9] [2]. The smallness of the mixing pa- 
rameter cosa is favored by the current LHC Higgs 
coupling measurements, and we will study the issue 
in later part. In this scenario, the coupling of the 
SM-like Higgs to fermions ypşf approaches to the 
SM one, thus the flavor-violating processes mediated 
by the SM-like Higgs boson are almost completely 
suppressed. 


II. SEVERAL RELEVANT TOPICS 
WHICH ARE SPECIFICALLY ADDRESSED 


A. Constraints on the parameter space of the 
afore model 


At first, we explore possible constraints coming 
from B physics and Electroweak precision test, and 
find that the model is more advantageous over the 
Type-II 2HDM as it may evade those constraints 
in the situation of my+ ~ ma because the tan 8 
enhancement effect does not exist. 


1. Constraints from B Physics 


A light a can mediate the initial and final states 
of decay Bs > u” uT in addition to the SM contri- 
bution, hence imposes a stringent constraint on the 
model. For Ma < mz, the correction due to an a 


exchange at the s-channel was calculated and 
the results are 
BR (Bs > utp") ~ BR (B, > HU) om (17) 


2 
vmpB, s2 p”? phe f (Tt, Yt, r) 


4mlm, M2) Y (a) 


xil + 


for œ? = 0, and 


uU 


BR (B; > wtp) ~ BR (Bs > utp), (18) 


2 
vmpB, s2 pt ph f (Xe, Yt, r) 


14 
my(my, =Ma) Y (2) 


x 


for pt = psg, with z; = me ime, Yt = ie | Wha 
r = m3,./miy, and the f and Y functions can be 
found in Eq. |A.3] of the Appendix. 

The average of the LHCb and CMS measurements 
on this mode is BR(B, > uty) = (2.9+0.7) x 
107° 54|. This could be compared against the 
SM prediction, which is taken to be (3.65 + 0.23) x 
107° [56]. We note that the tan 8 enhancement 
effect in the calculation of Bs > utu with the 
Type-II 2HDM [48] [51] [57] does not appear in our 
model and thus the constraint from B, leptonic de- 
cay is relaxed. 


2. The T parameter 


In the Type-III 2HDM, the constraint on the T 
parameter imposed by the electroweak precision test 
may be the most stringent. Following the method 


4 


proposed by the authors of Ref.{12} [58], the T pa- 
rameter in our model is obtained as, 


_ 1 2 
= 2 pa 
lrs miy 


x (F(m2,.,,m3,) — cos? OF (m3,m3))| a9) 


[F (mi, m+) cos? 0 + sina 


with 


The a/A and h/H mixing is highly constrained 
by the current LHC data to be around sina ~ 1 
and cos@ ~ 1, along with the parameter choice 
my+ ~ ma, the T parameter is suppressed, as seen 
in Eq. therefore it does not actually affect appli- 
cations of this model as indicated in [59]. 


B. Relevant processes under investigation 


In the following parts, we investigate several rel- 
evant processes in this extended Type-III 2HDM 
model. All the puzzles about the h > pT excess, 
muon g — 2 discrepancy, constraints coming from 
T — py, dark matter relic abundance and GCE ex- 
planations which were not solved in previous Type- 
HI 2HDM will be addressed. 


1. h — uT excess 


Existence of the flavor-violating Yukawa coupling 
in the extended Type-III 2HDM may possibly ex- 
plain the h —> uT excess observed by the CMS col- 
laboration. Now let us compute the branching ratio 
of h > pr in terms of our model, the result is shown 
as 


Mh 


(lurral? F lynurl?) ’ (21) 
where I» is the total decay width of the SM-like 
Higgs boson. To meet the observed excess, the flavor 
mixing should be of a magnitude 


BT \2 ou J2 /T 
pt =, / 10 tb | ~ ooo (4). (22) 


For T'a, ~ 4.2 MeV and cga = 0.01, as long as 
p'" ~ O(0.1) is reached, the new model is able to 
accommodate the h > uT excess. 


2. muon (g — 2), anomaly and T > uy 


The previous study indicated that new physics 
may contribute to the anomalous magnetic moment 
of the muon and radiative process T > py via a 
chirality flipping dipole operator [60 


ij 
(Heo Prej Fap +h.c., (23) 
NP 

where i,j denote the flavors of the external lep- 
tons, Fag is the electromagnetic strength tensor 
and o°? = i [v7 77]. The diagonal component 
contributes to the anomalous magnetic moment of 
muon, whereas the off-diagonal component corre- 
sponds to a dipole transition from T to u. It is 
noted that the coefficient £ is derived in var- 


ious new physics models and has different values 
which would help to determine the corresponding 
parameter space. In this work, we are going to de- 
rive this coefficient in the extended Type-HI 2HDM. 
The flavor-violating Yukawa couplings and the newly 
introduced pseudoscalar induce additional contribu- 
tions to (g — 2),, and T + py via one-loop and two 
loop Barr-Zee diagrams [61], as shown in the two 
panels of Figure. We include these extra contri- 
butions into the numerical computations. 


1. The muon g-2 anomaly 
In our case, the model-dependent coefficient 
C*i / A3, p could be expressed as 

Re{C"} vw 
Ayp V2 


ebay, (24) 


Amy, 7 


The contributions coming from the two CP- 
odd Higgs bosons give rise to 


ar, loop ~ de Sar ln ms 3 (25) 
Yor Srm m? 2 


a f ee 


(26) 


wee 


dal = 
ah An3m 


where ¢ = ve a, rf = m4 /m⁄, and yor fr is de- 
fined in Eq. (16 TO), fe could be found in Eq. = 
of the attached Pa. The W and Go. 
stone loops would not contribute to the eed 
Zee diagrams for Ca(ayww = 9. 


2. The T + py process 


With the lepton flavor violating Lagrangian 
L=emAgl; lio!’ q (APP, + AŻ Pr); (27) 
+h.c, 
the coefficients C* in Eq. (23) could be ex- 


pressed in terms of the form factors Az and 
Ar as done in Ref. [60], 


CT = em, An pAR Cut Z em, An pAr" 
Vw V2v 
(28) 
Thus, the branching ratio of T > puy is calcu- 
lated as 
BR(t > py) = BR(T > pw) x (29) 
48 em 
Aor Gem (Ar|? + |Aal?)- 
Gh 
$ 


t/u t/u 4 


FIG. 1: The Feynman diagrams for r > uy and (g—2),, 
where ¢ represents h, H, A,a and the loop at the bottom 
panel corresponds to t, W and Goldstone loops. 


Due to |Yhru| = |Ynur| in the model, we have 
|Ar| = |Ar| = |A|. The two CP-odd Higgs 
bosons contribute to the form factor A through 
1-loop and 2-loop Barr-Zee diagrams as shown 
in Fig. [I| and the new form factor is obtained 


as 
t,b 
A= = (Ar + Ad 2 (30) 
where 
2 
z Your Yorr ms 3 
Ai = “> T (az a (31) 


N.Q} Oem fo (rf) 


T Mmf 


Ay’ = 2Y YourYoff - (32) 


of 


The other contributions from the CP-even 
Higgs bosons are induced by the 2-loop Barr- 
Zee diagrams where an intermediate photon 
and a W-boson are involved [12] [62]. 


3. DM annihilation and GCE 


For Ma K ma, the Dirac DM fermions annihilate 
into bb primarily through exchanging a at s-channel, 
the annihilation cross section for the relative velocity 
Up is given as 

pie? Ke (0 7)? som 
16m (s — m2)? + m2T2’ 


(33) 


where ys is the center-of-mass energy of the anni- 
hilating DM fermions, Mma and Ta are the mass and 
decay width of the mediator boson a respectively. In 
the non-relativistic approximation s ~ 4m2 + m%v%, 
thus Eq. [33] can be rewritten as 


p 3 TAAST 
~ 256am2 (6 + 2/4)? +72 


(ovr) (34) 


where y = Maľl'a/ 4m? and ô are two dimensionless 


parameters, and the kinematic factor ô is defined as 
ô = 1 — m? / (4m2). 

As long as 6 is not too small, the DM annihi- 
lation could occur in the region far away from the 
resonance, then the cross section is almost indepen- 
dent of the velocity. In this case, the GCE and cor- 
rect thermal DM relic density could be accommo- 
dated simultaneously provided that the parameter 
y2 (7)? 859 is adjusted to an appropriate value. For 
small 6, the resonance effect would enhance the DM 
annihilation cross section. For that case, adjusting 
the parameter y% (p7 )?859/(d° + 77) can give a rea- 
sonable explanation to the GCE observation. When 
ô > 0, Eq. [34] indicates that the magnitude of (ovp) 
decreases as the temperature increases and the pro- 
cess Xx — a — bb does not sufficiently reduce the 
DM abundance at the freeze-out epoch, therefore 
some other annihilation processes which affect the 
DM relic abundance must exist in the Lee-Weinberg 
evolution equation. 


IV. A SYNTHESIS OF ALL THE 
INGREDIENTS 


In this section, we perform a numerical analy- 
sis to investigate the CMS h — wr excess, muon 


(g — 2), anomaly and the T > py, as well as DM 
relic abundance and GCE in the extended Type- 
II 2HDM. The model has been implemented in the 
program FeynRules|63], and the model file of the 
form CalcHEP [64] has been employed in the pack- 
ages micrOMEGAs 4.1.8 [65] to calculate the relic 
density and the annihilation cross section of DM. 


The Higgs masses in the model are set as: my = 
150 GeV, my+ = ma = 300 GeV. Here the value of 
myy+ is allowed by the flavor physics constraint [66], 
and m4 = my+ is suggested by the T parameter 
constraint as indicated in section. and the 
value of my employed in our computations is con- 
sistent with the vacuum stability requirement [67]. 
The magnitude of the two CP-even Higgs mixing 
angle a, is severely constrained by the recent Higgs 
data [68], and its closeness to 7/2 will be explored 
in later parts of this work. 


The invisible and undetected decays of SM-like 
Higgs boson are accounted as the decays of beyond 
SM (BSM) Higgs boson [68], and it is denoted as 
T(h > aa) +T(h > wt) = Tpsm(h). Since the 
contribution of BSM to the decay branching ratio is 
bounded bellow 0.34 at 68% CL, the decay of h > aa 
and h — wuz in this model would be constrained. 
The best fit of the branching ratio of the undetected 
decays is < 0.23 at 68% CL, and this limit allows 
that of BR(h > pr) = (0.8479:32)%, therefore it 
implies that the constraint on the flavor violating 
process is relaxed. As Ma < mp,/2, the a/A mixing 
angle 0 dominates the exotic decay rate of the SM- 
like Higgs h > aa and changes the total decay width 
of the SM-like Higgs, by which the prediction value 
of BR(h — wr) in this model would be affected. 
The magnitude of 0 is required to be at order of 
O(0.01 — 0.1) required by the present Higgs signal 
fit which is also welcome by the estimate of the DM 
relic abundance and GCE interpretation since this 
value affects the magnitude of (cv,)( see Eq34). 

The parameter dependence of the coupling rea 
denoted as pə for simplicity) needs to be carefully 
analyzed for the following reasons: the coupling 
Ppp is responsible for one of the dominant decay 
modes of the SM-like Higgs boson h —> bb, thus 
any change of py, would affect the theoretical pre- 
diction on BR(h > uT); meanwhile the DM annihi- 
lation xx — bb is supposed to be the dominant one 
and its cross section is proportional to the square of 
Pbpb, aS given in EqB4] Pop is set at the same order 
aS pur, ie. as aforementioned in sec|IIIB1} The 
contributions of tau lepton and top quark dominate 


BR(t > uy), as shown in Eq]31] and Eq[32] Thus 
the two parameters p77 and ptt (denoted as p+, and 


ptt for simplicity) need to be explored, here we use 
the parameter range of p;-(pit) ~ Poo. Even though 
we do not take the Cheng-Sher ansatz for p} we still 
choose a negative value for p? as in the Cheng-Sher 
ansatz [9] 02], except for p+- and p while consider- 
ing the current experimental constraint of T > py. 

To interpret GCE, the DM fermion mass and the 
coupling of the DM fermion to the mediator (a, A) 
are set as my = 30 GeV and yx = 0.5, and the range 
of Ma ranges between 30 GeV-95 GeV to account for 
the resonance effects in the DM annihilation process. 

To obtain the parameter spaces favored by the 
physical picture including all the aforementioned 
constraints, we perform a complete numerical anal- 
ysis for all possible parameter spaces: a— pur, &— 0, 
0 — Pur, Pob — Pur, Ptt — Prr, and Ma — 0, with rele- 
vant parameters being free within the ranges of: 30 
GeV< Ma < 95 GeV, 0.025 < 0 < 0.1, 1.475 < 
a < 1.57, —0.115 < pur < 0, —0.3 < pæ» < 0, 
—0.28 < pu < 0.28, —0.05 < prr < 0.05 based on 
aforementioned arguments. For each specific param- 
eter space listed above, several parameters need to 
be fixed as shown in the following 


e for parameters spaces of a — 0 (œ — pur): 


Pup = —0.01, pr- = 0.012, pr = —0.2, Pob = 
—0.2, and Mma = 50GeV(46 GeV), pur = 
—0.102(0 = 0.06) ; 


e for parameters spaces of pit — Prr (Pob — Pur): 
Ma = 50GeV, a = 1.546, puu = —0.01, 0 = 
0.045, pop = —0.2 and pur = —0.102 (prr = 
0.012 and pu = —0.2); 


e for the pur — 0 (Ma — 0) parameter spaces: 


a = 1.546, pr- = 0.015, pu = —0.2, py = 
—0.2, and pu = —0.01(—0.02), ma = 50GeV 
(pur = —0.102). 


Conducting a numerical analysis by means of the 
above parameter setup, the relevant processes are 
depicted by the figure | (see the caption of the figure 
for details). 


1. The h > pr 


To understand the CMS excess BR(h > uT) = 
0.847): 32%, the relation between the coupling 
Pur and a given in Eq. [22|should be satisfied. 
The region colored by magenta in the top-left 
plot of Fig. [2] is allowed to explain this ex- 
cess. With a properly fixed a, the both plots of 
the middle panel show that a sizeable coupling 
Pur ~ 0.1 is required. 


The exotic decay mode h — aa with a large 
branching ratio can efficiently change the to- 
tal Higgs decay width, thus the a/A mixing 
angle 0 and the mediator mass ma affect the 
prediction of BR(h — ur), as plotted in the 
bottom-right panel. 


The whole magenta region in the bottom-right 
panel with the value of BR(h —> uT) = 0.61 
implies that when @ is relatively small or in a 
case of 2m, < Ma, the decay process of h + aa 
cannot occur, so would not affect the Higgs 
total width and becomes irrelevant to h > uT 
process. Eq[2]] clearly interprets the situation. 


It is worth indicating that when this work is 
close to be finalized, the ATLAS collaboration 
published new analysis on h > uT, they ob- 
tained a slightly smaller excess compared to 
the CMS result, while its upper bound is con- 
sistent with the CMS result [69]. 


. The muon g — 2 anomaly 


The contribution of two-loop Barr-Zee dia- 
grams to ĝa, is negligible because of the small- 
ness of cosa. The dominant contributions to 
ôa, include the one-loop diagrams where CP- 
even Higgs H and CP-odd A are mediators, 
especially, the two diagrams respectively pro- 
vide negative and positive contributions. 


The contribution of the one-loop diagram 
where CP-odd a is the mediator to ĝa, is pos- 
itive, it becomes larger for smaller Mma and/or 
larger 0, therefore cancels out the contribution 
of the diagram where H is involved. The situ- 
ation is depicted at the top- and bottom-right 
panels of Fig. |2| (by the red dashed curves). 


To gain the outcome result which is con- 
sistent with the present experimental value 
(ôa, /107° = (2.61 + 0.8) [70]), a sizeable p,- 
is required, it is also consistent with the CMS 
h — pr excess, and its dependence on pur 
is demonstrated in the plots(top-left and the 
middle panels). 


Within the allowed parameter space in p+r- 
prt, a value of muon g-2: da, = 2.80 x 107° is 
reached. 


. Constraints from T —> uy 


The dominant contributions come from the 
one-loop diagrams (A-loop and H-loop di- 
agrams) and the Barr-Zee diagram with 
top quark being involved, as shown in 
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FIG. 2: The corresponding results are presented as functions of a, pop, Ptt, Prr, Pur, Ma and 0. The cyan (orange) 
regions give the DM annihilation cross section of (ovr) = (0.5 — 4) x 107*%em?/s((0.5 — 1) x 107cm? /s). DM relic 
density Qh? (= 0.1197 +0.0022) is depicted by blue lines and contours. Magenta region/contours are for CMS h > uT 
excess. Red dashed lines are contours of ĝa, /107°. The black curves are contours of BR(t —> p)/107%, and the red 
region in the bottom-left panel is consistent with the experimental constraint < 4.4 x 1078. The green curves are 
contours of BR(h > aa). The gray region is ruled out by the constraints of CMS Higgs signal strengths. 


Eq. (30]31]32). Existence of opposite signs be- tween the contributions of the CP-even Higgs 


H (negative) and the CP-odd Higgs A (posi- 
tive) one-loop diagrams [12] leads to a cancela- 
tion effect, however, it does not occur as long 
as ma > my. With a large coupling pu and 
under the limit sina ~ 1, the Barr-Zee dia- 
gram involving top-quark composes the domi- 
nant contribution (positive) to T > py. 


Additionally, like the case of ĝa, as Ma is 
smaller and @ is larger, the importance of the 
one-loop where a is involved to BR(r > uy) 
enhances as shown in the bottom-right panel 
of Fig. As well, a smaller œ would also in- 
duce an enhancement of the contribution of the 
h mediated Barr-Zee diagram with a W loop. 
Thus due to the contributions of this Barr-Zee 
diagram and the a one-loop diagram, for a re- 
gion with larger 0 and smaller a, our calcula- 
tion would predict an even larger BR(t > uy) 
(see the top-right panel). 


The theoretical prediction on BR(t > puy) 
with respect to p,,-(plotted at the top-left and 
middle panels) shows that its behavior is sim- 
ilar to that for da,. 


It is noteworthy that since the the form fac- 
tor A is related to tau lepton and top quark 
through the concerned loops, (see Eq. (31]32)), 
the current experimental bound BR(t —> 
uy) < 4.4 x 1078 [ZI] [72] constrains p,- and 
pit Strictly. The red region in the bottom-left 
panel of Fig. |2| is allowed. It is interesting 
to note that pu and p+- should have opposite 
signs as favored by the data. This result agrees 
with that given by the authors of Ref. [II]. 
One can see that the upper bound demands 
pr- to be small as about |p,-| < 0.04. 


. DM relic density and GCE 


As long as assuming that the observed GCE is 
caused by dark matter annihilation, the lower 
bound of the annihilation cross section (ovp) 
should be about 0.5 x 10~76cm3/s as discussed 
in Ref.[22]. The upper bound of (cv,.) is deter- 
mined to be 4.0 x 107° cm?/s at 95% CL. [33]. 
Also the data of Pass 8 of Fermi-LAT from 
dwarf spheroidal satellite galaxies set a new 
upper bound (~ 1.0 10~76cm3/s) on the dark 
matter annihilation cross section at the dark 
matter mass 30 GeV. The areas in Fig. [2] that 
give an annihilation cross section of (av,;) = 
0.5 — 4.0(1.0) x 107?ĉcm?/s are depicted as 
the cyan (orange) regions. As for the top-left 
(bottom-left) parameter space, the calculated 


values of (ovp) = 1.8(2.5) x 107cm? /s. 


The blue curves and contours from Fig. [2|rep- 
resent the correct dark matter relic density 
(Qh? = 0.1197 + 0.0022 [74]). We note that 
the range of (cv,.) favored by the DM relic den- 
sity and GCE as expected highly depends on 
Ma, 9, and po (see Sec. [IIB 3). Due to the 
enhancement of the annihilation cross section 
when the mediator mass m, is close to ~ 2m,, 
the dark matter relic density rules out a range 
in the parameter space (see the bottom-right 
panel). 


With a sizeable pẹ ~ O(0.1), the DM relic 
density and GCE could be tolerated simul- 
taneously in a range of (ovr) = (0.5 — 4) x 
10- “em? /s, 

However, the results newly reported by the 
Fermi-LAT and DES collaborations constrain 
(our) to be smaller than 1.0 x 10~?°cm3/s 
for the 30 GeV dark matter fermions, which 
is slightly smaller than the value required by 
the thermal relic abundance, thus there should 
exist other additional DM annihilation chan- 
nels to make up the correct DM relic den- 


sity [75] [76]. 


. The CMS constraints 


Recently, the CMS collaboration combined 
the comprehensive sets of production and 
decay measurements for the 125 GeV- 
Higgs boson, including decay channels into 
yy, Z2Z2*,WW*,rt+r-, bb, and utu, and 
found no significant deviation from the stan- 
dard model predictions [68]. That synthesis 
should severely constrain the parameter spaces 
of all built models. 


More specifically, the signal strengths for 
T*7~, bb, WW*, ZZ*, yy channels are defined 
as 


a a(h) x BR; 
Ki GSM (h) x BRY ’ 


which can be gained by fitting the CMS data 
within lo tolerance [77]. In the text o(h) and 
o°M(h) ( BR; and BRIM) correspond to the 
Higgs production cross section (decay branch- 
ing fractions of the five decay modes) pre- 
dicted by respectively this model and the SM 
as BR; = T(h =. it) /(T gsm (h) + T'gm(h)). 

We find that the upper bound on pw and 
the lower bound on py set more rigorous con- 
straints on the parameters of this model, and 


(35) 


10 


CMS[63]_ | ma-0 | a-0 | a— pyr | pur—O | prr- pit | Pur — poo 
u. | 0.91 £0.28 | 0.94-1.07| 0.92-1.11 | 0.93-1.10 | 0.94-1.08 | 0.85-1.19 | 0.92-1.0 
i | 0.84 £0.44 | 0.58-0.67| 0.63-0.65 | 0.64-0.65 | 0.58-0.67 | 0.64-0.69 | 0.64-0.66 
uwl| 0.83£0.21 | 0.88-1.0 | 0.88-1.04 | 0.89-1.03 | 0.88-1.0 | 0.96-1.04 | 0.87-1.03 
juz} 1.0£0.29 | 0.88-1.0 | 0.88-1.04 | 0.89-1.03 | 0.88-1.0 | 0.96-1.04 | 0.87-1.03 
jy | 1.12 £0.24 | 0.88-1.0 | 0.88-1.04 | 0.89-1.04 | 0.88-1.01 | 0.97-1.04 | 0.88-1.04 


the grey areas in Fig] are excluded. With the 
allowed spaces of the parameters, the ranges 
of the signal strengths u; are given in Table. 


The CMS Higgs decay signal strength strongly 
constrains the h/H, a/A mixing angles a and 
0, concretely it demands both sin a and cos 0 to 
be close to unity (see the middle-left and top 
panels). The constraint on a is found to be 
compatible with the parameter space obtained 
by accounting for the h —> uT excess. ppp is 
also severely bounded since it could affect the 
primary decay mode of the SM-like Higgs sig- 
nificantly. As is shown in middle-right panel 
of Fig. |2| which is a contour diagram of pbb- 
Pur, One notes that only a small band with pr» 
being between —0.15 and —0.22 is not ruled 
out. 


Meanwhile from the contour diagram of ma—0, 
one can see that as long as BR(h — aa) is 
larger than 0.12, the grey regions in the fig- 
ure are excluded by the CMS Higgs signal 
strengths. 


For the not-yet-detected channel h — aa , 
the pseudoscalar primarily decays into bb, then 
there should be an additional contribution of 
the process h — aa — 4b in the h — bb 
searches [78], and the data of the CMS exper- 
iment would definitely constrain the 
coupling between hbb as long as Ma < mp/2. 
The upper bound on the undetected decay of 
Higgs will be further improved as the bb pair 
production is measured at the 13 (14)TeV. 
Moreover the h — 2b2v searches at LHC [48] 
could give more rigorous constraints on the pa- 
rameters Ma and ð. 


. Collider observation of the pseudoscalar 


For the benchmark scenarios of the work, when 
Ma > 2my,, the masses relation mA > Mp +Ma 
opens the channel pp(gg) > A > h(h > 
yy)a(a — xX) for the observation of the pseu- 


V. 


TABLE I: A comparison of the signal strengths in the allowed parameter spaces of this model with the Higgs data 
coming from CMS [68]. 


doscalar a at the LHC, which is found to be the 
case of the mono-Higgs searches in [80]. The 
backgrounds are dominated by the SM pro- 
cess pp > Zyy! with Z > vv, and the Higgs 
associated production process pp > Zh with 
Z — vv. The collider analysis is performed by 
generating signal and background events with 
MaApGRAPH5_AMC@NLO [81] [82] at 14 TeV, 
and then passing on to PYTHIA 8.1[84] for par- 
ton shower and hadronization, and the detec- 
tor simulation is conducted by DELPHES 3[85] 
at last. After implement the selection cuts 
Myy € [120,130] GeV, as well as Hr, PX > 76 
GeV following Ref. [80], the number of events of 
signal and background are obtained as 37 and 
48 respectively, thus the significance for the 
observation of a is found to S/V S + B ~ 4 for 
the benchmark of ma — 0 with 0 = 0.08, mg = 
76 GeV and with an integrated luminosity 
L = 300 fb7}. 


The probe of the pseudoscalar could also be 
conducted with hard b-jets and large missing 
transverse energy [86] when ma > 2m, and the 
pseudoscalar dominantly decay to dark mat- 
ter. And when the pseudoscalar decay mostly 
to b-quarks, the collider search of the pseu- 
doscalar could be found in Ref. [24]. 


CONCLUSIONS AND DISCUSSIONS 


The discovered anomalies by the LHC experi- 
ments and observations of the Dark matter greatly 
excite the curiosity of human beings and inspire en- 
thusiasm of searching for new physics beyond stan- 


dard 


model. However, so far, the trend is not 


1 Here, we would like to mention that including the back- 
ground Zy + jets (with a jet faked a photon) might reduce 
the observation significance a bit 


very successful. Even though we know new physics 
must be around, but do not know its scale. Direct 
and indirect search for dark matter, LHC experi- 
ments, long-baseline and short-baseline neutrino ex- 
periments and numerous lower energy experiments 
including BELLE, BES and many others provide 
hints towards new physics beyond standard model, 
however at the same time, set more and more rig- 
orous constraints on those models which have been 
proposed to explain the anomalies observed in as- 
tronomy and high energy experiments. Some of the 
models survive the so-far measurements and many 
have been ruled out. 

In this work, we have extended the Type-III 
2HDM by introducing a pseudoscalar a which can 
mediate the DM pair annihilation process. In a re- 
cent paper, Han et al.[87] also extended the 2HDM 
with an aligned Yukawa sector to explain the (g—2),, 
excess, in comparison, our scheme is somewhat dif- 
ferent from theirs. We not only consider the (g—2),, 
excess, but also many other constraints from both 
earth experiments and cosmology. In this frame- 
work, the LFV process h —> pt observed at LHC 
is addressed, the possibility of explaining the muon 
g—2 anomaly, dark matter relic abundance and GCE 
have also been investigated, the role played by the 
newly introduced a is studied in some details. It 
is found that there indeed exist certain parameter 
spaces which can tolerate h + uT excess, muon g—2 
discrepancy and dark matter relic abundance. 

With the flavor violating coupling ppr ~ O(0.1) 
and a tiny mixing between h and H around sina ~ 
1, the observation of the h —> yr excess could be 
easily accommodated. The pseudoscalar a opens an 
important undetected decay channel for the SM-like 
Higgs (h > aa) as Ma < mp/2, thus affects the 
branching ratio of h > Tpu. It also plays a role to 
explain the discrepancy between theoretical predic- 
tion and data of (g—2),,. A smaller ma and a slightly 
larger CP-odd a/A mixing 0 helps to interpret the 
muon g — 2 anomaly. Increasing the CP-odd Higgs 
mixing angle and |p,,-| increases BR(t > uy), thus 
would further constrain the parameters space of our 
model. The measurement of the branching ratio of 
T — py sets a stringent bound on the flavor con- 
serving Yukawa couplings p and p+r, moreover, it 
determines an opposite sign between the two cou- 
plings. There are parameter regions in our model 
allowed by the current experimental data of T > uy 
where both the measured h — uT excess and the 
muon g-2 anomaly can be explained. 

To account for the dark matter relic density, the 
a/A mixing angle should be of order O(0.01) and 
Pop ~ O(0.1) is required, however this parameter 
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region does not coincide with that favored by the 
GCE observation if it is postulated that the GCE 
is fully coming from the yx DM pair annihilation. 
This inconsistency may imply that there are other 
sources to result in the GCE besides the pure xX 
DM annihilation mechanism. 

The LFV process h — yt is studied in terms of our 
model and our prediction on BR(h — pT) is quali- 
tatively consistent with that observed by CMS and 
ATLAS at 8 TeV, however to make a decisive conclu- 
sion, more data are needed and the LHC Run II of 14 
TeV should help. A synthesis of data accumulated 
by high energy collider LHC, the future SPPC of 
100 TeV and maybe some lower energy experiments 
as well as the new astronomical observation would 
make the whole picture clearer, then we will be able 
to judge whether this model indeed works or needs 
to be further modified. For the case of ma > 2m,, 
the mono-Higgs search provide one possible probe of 
the pseudoscalar. 


Appendix: Loop functions 


The function fọ and h using for the calculation of 
g — 2 and T > uy are given by, 


Oor 1 1 x(1— x) 
faal) = sf Hra r 


The functions f (z,y,r) and Inami-Lim function 
Y (x) being used in Eq. are 


(A.1) 


x r(a—1)-—2 _ xlogr 
Htewr=3| e-e- Ete- 
ylogy x logy 

u- pae AA 
ro- [eh = (A.3) 
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